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Abstract. Effects of dipole-dipolc interactions on the magnetic relaxation have been 
investigated for three Fe-C nanoparticle samples with volume concentrations of 0.06, 
5 and 17 vol%. While both the 5 and 17 vol% samples exhibit collective behavior due 
to dipolar interactions, only the 17 vol% sample displays critical behavior close to its 
transition temperature. The behaviour of the 5 vol% sample can be attributed to a 
mixture of collective and single particle dynamics. 
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1. Introduction 

The dynamics of systems of interacting ferromagnetic nanoparticles has been the focus 
of extensive research in recent years. A question of considerable controversy has 
been the existence of a phase transition to a low-temperature spin-glass-like phase 
H, 0, 0, H || |], 0, H ||. Recent studies on strongly interacting systems have 
reported a critical slowing down ||, [| and a critical divergence of the non-linear 
susceptibility M. In addition, non-equilibrium properties such as ageing, memory and 
rejuvenation phenomena have been observed in magnetic relaxation and low frequency 
AC susceptibility measurements at low temperatures 0, lOf- These observations indicate 



the existence of a low-temperature spin-glass-like state. However, some questions about 
the collective state of interacting nanoparticles are still not resolved. In a recent 



experimental study of an interacting maghemite nanoparticle system ||11|| , it was not 
possible to find evidence of a finite temperature transition to a spin-glass-like phase. 
Still, the particle system exhibited non-equilibrium dynamics in magnetic relaxation 
experiments, typical of spin glasses. The absence of a critical behaviour was in this 
study linked to mixing of collective and single particle relaxation effects. Mixing of this 
kind can for instance originate from a broad particle size distribution, or it may be 
observed in a system with a small volume fraction of particles and therefore comparably 
weak inter-particle interactions. 

In this work, extensive studies of the magnetic dynamics in a close to 
monodisperse nanoparticle system are presented. Three samples having different 
volume concentrations of nanoparticles were investigated; one dilute and nominally non- 
interacting sample and two more concentrated interacting samples. Both interacting 
samples exhibit collective behaviour due to dipolar interactions at low temperatures, 
but only the most interacting sample displays critical dynamics in the investigated time 
window and shows equilibrium dynamics quite similar to that of an atomic spin glass. 



2. Experimental 

The particles were prepared by thermal decomposition of Fe(CO)5 (20.0 ml) in a mixture 
of carrier liquid (50.0 ml cis-trans decalin) and surfactant (4.0 g oleic acid) by the method 
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described in This method leads to the formation of surfactant coated particles of 
the amorphous alloy Fe^Ci^. (x ps 0.2-0.3). The decomposition was carried out slowly 
by using low heating powers. Great care was taken at all times to avoid oxidation 
of the particles and a flow of Ar/H 2 (98/2%) gas through the reaction vessel was 
maintained throughout the preparation. After preparation, the decalin carrier liquid 
was evaporated at 140 °C under reduced pressure using a gas flow to drive out the 
vapour and a small amount of de-gassed xylene was added such that the volume fraction 
of particles in the resulting base ferrofluid was about 5 vol%. More dilute samples 
were obtained by further dilution of this ferrofluid by addition of xylene and more 
concentrated samples were obtained by evaporating the xylene at room temperature in 
vacuum. All sample handling after the preparation of the base ferrofluid was carried out 
in an argon glove box. The chemical state of the iron in the ferrofluids was checked by 
Mossbauer spectroscopy, which showed that no significant oxidation of the particles had 
taken place during the sample preparation and sample handling. The particle volume 
fractions in the resulting samples were estimated from the iron content determined 
by atomic absorption spectroscopy. The sizes of the particles were determined from 
transmission electron microscopy (TEM) studies. A droplet of a dilute ferrofluid was 
placed onto a grid, that was left in air for about a week to ensure full oxidation of the 
particles. The size determination from the resulting TEM micrographs was subsequently 
corrected for the change in density due to the oxidation of the particles. The shape of 
the particles was spherical to a good approximation. The particle size distribution was 
obtained using the computer analysis method described in |T3. The analysis of 1579 
particles yielded the average particle diameter d = 5.3 nm with the standard deviation 
0.3 nm. The volume-weighted volume distribution was well described by the lognormal 
distribution, f(V)dV = (2iry 1 / 2 (aV)- 1 exp[- \n 2 (V/V m )/2a 2 ]dV, with V m = 78.2 nm 3 
and a = 0.13. The saturation magnetisation was estimated to 1 x 10 6 Am -1 from 
Langevin fits to magnetisation vs. field curves at several temperatures well above the 
superparamagnetic blocking temperature. The present preparation batch is slightly 
different from that studied previously in @, [| We have studied three samples with 
particle concentrations of 17 ± 4, 5 ± 1 and 0.06 ± 0.02 vol%. In the most concentrated 
sample, most of the carrier liquid was evaporated and the sample had a paste-like 
consistency. This resulted in poor thermal contact between the sample and the sample 
container. Thermal stability was not achieved until several minutes after a change in 
temperature, preventing systematic studies of the non-equilibrium dynamics for this 
sample. 

Two experimental set-ups were used for the magnetic studies. A commercial 
susceptometer (LakeShore model 7225) was used for DC field scans and AC susceptibility 
measurements in a frequency range of 15 - 1000 Hz. A non-commercial superconducting 
quantum interference device (SQUID) magnetometer [nj was used for AC susceptibility 
measurements in a frequency range of 10 mHz - 9.1 kHz, zero field cooled (ZFC) 
relaxation measurements and magnetic noise measurements. The applied AC and DC 
fields were chosen small enough (0.1 Oe) to ensure linear response from the samples. The 
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background field was less than 1 mOe. The ZFC relaxation measurements were carried 
out by cooling the sample to the measuring temperature T in zero field, equilibrating 
the system for a wait time t w , and then applying a small DC field h and measuring 
the magnetisation M(T, t) as a function of time, t, after that the field was applied. 
For a slowly relaxing system, it can be shown that M(T, t) and the relaxation rate 
S(T,t) = h- l dM(T,t)/d\nt relate to the real and imaginary components of the AC 
susceptibility as 

M(T,t)/h^ X '(T,u), (1) 
S(T,t) «V(:Z», (2) 

7T 

with t = 1/uj. Magnetic noise measurements were performed in zero external field 
using an HP35670A dynamic signal analyser. The power spectrum of the magnetic 
fluctuations was measured in three overlapping frequency intervals (i) 0.003-12.5 Hz, 
(ii) 0.25-400 Hz and (iii) 8-12500 Hz. The background spectra were subtracted from the 
data. Data of the same order of magnitude as the background signal were not used in 
the analysis. The fluctuation-dissipation theorem ]jnj relates the noise power spectrum 
to the zero-field limit of the out-of-phase AC susceptibility as 

P(T,Lu) = 4k B T X " i - T,UJ \ (3) 

A comparison between the out-of-phase component obtained from AC susceptibility 
measurements and results obtained from zero field noise measurements verifies that 
equation (|3|) is obeyed and hence that the AC data are obtained in the linear response 
regime. 

3. Results and discussion 

3.1. General behaviour 

Figure [l] shows the AC susceptibility measured in the commercial set-up for the three 
samples. The most dilute sample with a 0.06 vol% fraction of particles is intended 
to serve as an experimental reference for a non-interacting system. However, the 
peak height of x" f° r this sample increases slightly with increasing frequency. This 
trend is opposite to the expected behaviour for non-interacting nanoparticles |18| and 



thus indicates that the dynamics of this sample is influenced by weak inter-particle 
interactions. Despite this, it is possible to assess the distribution of energy barriers 



with reasonable accuracy making use of the method described in ||19|| . This analysis 
yields r = 1 x 10~ 12 s and K = 0.9 x 10 5 Jm -3 . The extracted values of K and r 
compare reasonably well with previous estimates. For 5.0 nm Fe-C particles, K has been 
estimated to 1.3 x 10 5 Jm -3 , and it has been found that K increases with decreasing 
particle size to 3 x 10 5 Jm -3 for 3.2 nm particles |20, 21]. Studies of different preparation 



batches have reported values of tq in the range 2 x 10~ 12 s to 3 x 10~ n s 0, II, pi] , 

Figures ||] and ||| show the AC susceptibility data measured in the non-commercial 
SQUID for the 5 and 17 vol% samples, respectively. For the 5 vol% sample, ZFC 
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relaxation data are also included. A few differences in the behaviour of these samples 
as compared to the dilute sample can immediately be noticed. The peak height of x" 
increases significantly with increasing frequency, and the width of the peak increases with 
increasing frequency. The onset of a non-zero x"i is shifted towards higher temperatures 
and it becomes sharper with increasing volume fraction of particles. 

Figure [| shows the relaxation rate for the 5 vol% sample measured for two different 
wait times at different temperatures. Below 40 K, a clear wait time dependence is 
observed indicating that non-equilibrium phenomena plays a key role for the dynamics 
at low temperatures. The same sample has recently been subject to a detailed study 
of the low temperature non-equilibrium dynamics in which the memory effect in the 
ac susceptibility was observed and further characterised by temperature cycling ZFC 
relaxation experiments [|Kj. It was shown that the non-equilibrium dynamics in this 



sample is governed by long-range collective behaviour and that well-known concepts 
for spin glasses such as chaos with temperature and overlap length are necessary 
to describe the observations. This implies that the magnetic relaxation of the two 
concentrated samples must be analysed in terms of collective dynamics. A model for 
weakly interacting particles, such as that proposed in f23fl , where the effect of dipolar 
interaction on the relaxation time is accounted for by introducing the thermodynamic 
averages of the local dipolar field in a rigorous expression for the single particle relaxation 
time, is only valid at temperatures much higher than the freezing temperatures of the 
two concentrated samples. Below, we briefly introduce the concepts of dynamic scaling 
and discuss the observed dynamics for the two concentrated samples in terms of critical 
dynamics. 

3.2. Dynamic scaling 

A signature of a continuous magnetic phase transition is the divergence of the correlation 
length, £, when the phase transition temperature, T g , is approached from above as 
£/a = e~ u , where a is the average distance between interacting moments, e = T/T g — 1 is 
the reduced temperature and v is a critical exponent. According to conventional critical 
slowing down, the longest relaxation time due to correlated dynamics, r c , is related to 
the correlation length as r c oc (^/a) z , where z is the dynamic critical exponent. Hence, 
for T -> T+ 

r c = T*<r w , (4) 

where r* is a microscopic relaxation time. According to the dynamic scaling hypothesis 
24| , for T — * T+ and t/r* 3> 1 , the spin auto-correlation function can be written in the 



scaling form |2| 

q(t)=t-V™Q{t/T e ), (5) 

where (3 is a critical exponent and Q(x) is a scaling function. Using linear response 
theory it is possible from this relation to obtain the complex susceptibility and derive 
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the scaling relation p6| 

where u = 1/t and G(x) is a scaling function. The asymptotic behaviour is G(x) oc x y , 
with y = 1 and y = (3/zu for small and large values of x, respectively. Using 
the asymptotic behaviour of G(x) in the limit ujt c — > 0, the following relation holds 



26, 27 



l x " {T ; UJ } oc e — +/3 K r i-/V- (7) 

implying that the left hand side of equation ([7]) at each temperature reaches a frequency 
independent plateau. 

The meaning of r* is that it is the relaxation time of the individual magnetic entities 
in the system. For spin glasses, r* ~ 1CT 13 s and is the fluctuation time of an atomic 
moment. For nanoparticles, r* can be assigned to the superparamagnetic relaxation 
time of a single particle of average size. In a dense system, the dipolar interaction may 
modify this relaxation time compared that of isolated particles. However, as a first 
approximation, it is reasonable to assume that r* is close to the superparamagnetic 
relaxation time of an isolated particle, which in the relevant temperature range for our 
studies can be approximated by the Arrhenius-Neel expression. Below, we compare 
two approximations: i) r* = constant, and ii) r* = To exp(KV m /kBT). The first 
approximation has been used in previous work [|, |], [| ||, |28f , but it is only a good 
approximation if there is little variation of the single particle relaxation time in the 
temperature interval used for the analysis. 

In AC susceptibility experiments, the slowing down of the relaxation time r c 
(equation (^) can be derived from the temperatures where an onset of dissipation occurs 
(freezing temperatures) as a function of the observation time uj~ 1 . We have considered 
two criteria for the onset. In the first criterion, the freezing temperature is defined as 
the temperature at which x"(T,u) attains 15% of its maximum value. In the second 
criterion, the freezing temperature is defined from the relation x'(Tf,uj) = 0.98x e q(Tf). 
Figure ^(b) shows the freezing temperatures for the two concentrated samples obtained 
from x' data. The superparamagnetic blocking temperatures estimated from the peaks 
of the x" data for the dilute sample are included for comparison. 

3.2.1. 17 vol% sample First, we consider the approximation r* = constant. Using the 
out-of-phase susceptibility data, a dynamic scaling analysis according to equation (f|) 
results in T g = 49.5 ± 2 K, zu = 10.5 ± 2 and r* = 1(T 7 - 7±1 s @. This analysis is 



performed for reduced temperatures in the range 0.16 <e< 0.37 and corresponds to 4 
decades of observation times. A similar analysis using the in-phase susceptibility data 
results in T g = 50.5 ± 2 K, zv = 9.5 ± 2 and r* = io~ 8 - 3±1 s. The values of T g and 
zv from the two analyses are in good agreement, and the value of r* is, considering the 
temperature interval used in the analysis and the estimates of r and K given above, 
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within the limits implied by the Arrhenius-Neel expression. The derived values of zu 
also compare well to values found in previous work on nanoparticles: zv — 11 ± 3 in 
Q and zv = 10.5 ± 3 in ||, but less well to zv = 7.0 ± 0.3 found in reference |J. 
A good agreement is also found comparing with canonical 3D Ising and Heisenberg 
spin glasses, for which zv = 8 — 10 []2B, 27j| . An analysis according to the full scaling 
relation, equation (||), using T g = 49.5 ± 1.5 and zv = 10.5 ± 2.0 yields data collapse 
for (3 = 1.1 ± 0.2 (see reference ||28|| , figure 3). This analysis is based on all available 
x" data with temperatures corresponding to e > 0.01 (T > 50 K). The value of (3 is 
in good agreement with (3 = 1.2 ± 0.1 obtained from a full scaling analysis of the non- 
linear magnetic susceptibility on a similar sample ||, but larger than typical values of 
(3 = 0.5 — 0.8 reported for 3D Ising and Heisenberg spin glasses [^6], |29|, |3(|. The value of 
(3/zv » 0.11, extracted from the asymptotic behaviour of G(x) for large x, is consistent 
with the derived values of /3 and zv. 

Next, we consider the effect of using T* = r exp(i^l / m /fc s T). It should be noted 
that this introduces one extra parameters that can be varied in the analysis, and we 
have therefore explored the critical behaviour for a variety of choices of To and KV m to 
evaluate the robustness of the estimates of zv and (3. Critical slowing down analyses 
have been performed for various fixed values of Tq and T g and the values of zv and KV m 
have been estimated and used in the analysis according to the scaling relation, equation 
(^), to extract the value of (3. From these analyses, it is found that both the critical 
slowing down analysis and the full scaling relation can be fulfilled for T g = 50 ± 2 K 
with r = 1 x 10~ n s. The estimate of KV m /kB depends slightly on whether the x' 
or the x" data are used in the critical slowing down analysis, and the corresponding 
parameter intervals are KV m /kB = 500 ± 100 K using x' (shown in figure ||) and 
KV m /k B = 650 ± 100 K using x" ■ F° r both analyses, the extracted critical exponents 
attain the values zv = 8.5 ± 2 and f3 = 0.9 ±0.2. Other values of r of the order of 
10 -11 s yield slightly different values of KV m /kB but the same values of the critical 
exponents. The values of KV m /kB and r are of the same size as the estimates for the 
dilute sample, KV m /kB ~ 510 K and r ~ 1 x 10~ 12 s, but it is stressed that an exact 
correspondence is not expected. The values of zv and j3 are slightly smaller than those 
obtained from the analysis assuming a constant r*. The decrease of zv is due to the 
fact that the temperature dependence of r c / r* is weaker when r* is allowed to vary with 
temperature. The reduced value of zv and the extracted smaller value of (3 leads to 
(3/zv fa 0.11, which is the same value as found from the analysis using a constant value 
of r*. The data collapse of e _/3 x"(T)/x e(? (T) vs. ojt c to a single function G(x) according 
to equations (f|) and (Q) using r = 10~ n s and KV m /kB = 570 K (the average of the 
values obtained using x' an d x" m t ne critical slowing down analyses) is shown in figure 
[| The scaling is of the same quality as that assuming a constant r*, and the asymptotic 
behaviours are the same (in agreement with the estimates of (3 and zv). 

Figure |7|(b) shows x" i^) / u Xeq as a function of uj for different temperatures T > T g 
in a log-log plot for the 17 vol% sample. The prediction of plateaus at low frequencies 
from equation (0) holds regardless of a temperature dependence of t„ and figure 0(b) 
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shows that the prediction is well confirmed for this sample. In the same figure, 
the behaviour calculated from equation (|6|) is also included using the experimentally 
determined scaling function, G(x), and extrapolating using its asymptotic behaviour. 

3.2.2. 5 vol% sample For the 5 vol% sample, an apparent scaling according to equation 
(H) with a constant value of r* can be obtained for temperatures corresponding to 
0.25 <e< 0.8 (6 decades of observation times) with T g = 35.1 ± 2 K, zv = 10.8 ± 1 
and r* = 10" 4 - 5±a5 s using x" data, and with T g = 34.6 ± 2 K, zv = 10.8 ± 1 and 
r* = io _4 - 7±0 ' 5 s using x' an d ZFC relaxation data (shown in figure [|). However, 
deviations from scaling are found at lower temperatures (e < 0.25). Moreover, it is 
not possible to obtain data collapse according to the full scaling relation equation @, 
with temperatures corresponding to 0.02<e<0.8 (or to 0.25 <e< 0.8), for any choice 
of T g , zv and f3. Based on the estimates of KV m and r for the dilute sample and the 
temperature range used in the analysis above, it is expected that r* ~ 10~ 9 — 10~ 7 
s. However, the derived value (r* ~ 10~ 46 ) deviate by orders of magnitude from this 
range. Furthermore, even introducing an Arrhenius-Neel temperature dependence of 
T*, it is not possible to obtain unambiguous parameters from the critical slowing down 
analysis and to fulfil scaling according to equation (||). A simple manifestation of critical 
dynamics, which is independent of a temperature dependence of r*, is that x' / 'oJXeq 
should settle on plateaus for small values of uj. Conferring figure 0(a), it is seen that 
this prediction is not fulfilled for the 5 vol% sample (except for T > 65 K), and this is 
a further indication of non-critical dynamics in this sample. 

Several factors may contribute to deviations from critical behaviour and a possible 
non- divergence of the correlation length. First, there may be physical clustering of 
the particles {i.e., regions with a higher density of particles than average). The 
effect of clustering is two-fold: At comparably high temperatures, where short-range 
correlations are relevant, the stronger inter-particle interaction in particle clusters will 
enhance the local correlation length and shift the dynamics to longer time scales than 
expected from the volume fraction of particles and a homogeneous particle dispersion. 
In addition, the small scale heterogeneity will create a dispersion of length scales for 
the collective dynamics that modifies and possibly limits the growth of correlations 
at lower temperatures. Second, the polydispersivity of the particle system leads to a 
distribution of single particle relaxation times with a width that increases significantly 
with decreasing temperature (and is thus more important for the 5 vol% sample than 
for the 17 vol% sample). At lower temperatures, the largest particles may therefore 
become thermally blocked on time scales comparable to the longest time scale related 
to the collective dynamics and act as random magnets instead of taking part in the 
collective dynamics. This may completely mask or even obstruct long-range collective 
behaviour, as discussed in [11]. The non-critical behaviour of the present 5 vol% sample 



is most likely both due to the formation of clusters of particles (see section below) and 
the dispersion of single particle relaxation times. 
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3. 3. Relaxation function 

For a spin glass and T <T g , S(T, t) is history dependent and has a non-trivial variation 
with t. However, when quasi-equilibrium is probed, one can write S(T,t) oc t~ v<yT \ 
where y(T) attains a positive value close to zero f25j. For T > T g , the time dependence 



of S(T, t) is determined by the critical dynamics. As discussed previously, the critical 
dynamics result in fluctuations spanning all time scales between r* and r c , and this 
gives a very broad spectrum of relaxation times. From equation @ and using the 
asymptotic behaviours of G(x) in equation (||), it is easily seen that S(T,t) oc t~ y ^ 
with y(T) = p/zv « 0.1 for < t < r c and y(T) = 1 for t > r c ]3l|. At intermediate 
observation times, y(T) monotonically increases with t between the two extremes and 
a cross-over (or 'knee') from an essentially flat S(T,t) vs. t curve to S(T,t) oc t~ l is 
expected at t ~ r c . 

Figure [| shows the relaxation rate vs time for the three samples at different 
temperatures in the range 20 K < T < 70 K. For the 0.06 and 17 vol% samples x" data 
are shown, and for the 5 vol% sample results from ZFC relaxation, AC susceptibility and 
magnetic noise measurements are included. For the 17 vol% sample, x" data for T > T g 
obtained from equation (^j) using the experimentally estimated G(x) and the expected 
asymptotic behaviour for small x are also shown. In the low-temperature region, T < 30 
K, the magnitudes of the relaxation rates for the two interacting samples are much 
smaller than that for the 0.06 vol% sample. This is a signature of collective dynamics 



il| , p2| . Moreover, the relaxation rate of the 17 vol% sample is smaller than that of 
the 5 vol% sample. For temperatures T > 40 K and in the investigated time window, 
the relaxation rate of the 0.06 vol% sample decays with increasing observation time, 
while the relaxation rate of the the 5 vol% sample at short time scales exhibits a weak 
frequency dependence followed by a knee and a decrease towards zero at longer time 
scales. This approach towards zero relaxation rate (y(T) ~ 0.5 — 0.6 between 45 and 
55 K) is slower than that of an atomic spin glass (for which y(T) = 1) and consistent 
with the lack of plateaus in the x" '0^0 / '^Xeg curves for the corresponding temperatures 
in figure 0(a). The time dependence of the relaxation rate at time scales shorter than 
this knee is also uncharacteristic of a spin glass as it is non-monotonic in time and 
shows a broad maximum. For the 17 vol% sample and in the temperature range 35 K 
<T<50 K, the relaxation rate exhibits a slow, monotonous decrease with increasing 
observation time. The relaxation rate follows a S(T,t) oc t^ y ^ dependence, with y(T) 
decreasing with decreasing temperature, mimicking the expected behaviour of a spin 



glass system |2"5] . It should be noted though, that the particle size distribution will 
cause a temperature dependent and broad onset of the response on short observation 
times, as can be envisaged by the relaxation rate curves at the lowest temperatures. 
At temperatures above T g , T > 50 K, a knee appears in the relaxation rate also for 
this sample and the approach towards zero follows a S(T,t) oc t~~ l . At temperatures 
T > 65 K, the cut-off in the relaxation rate is equally sharp for both interacting samples, 
and the dynamics of the two samples are rather similar. This indicates, for the 5 vol% 
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sample, that the particle moments participating in the dynamics on long time scales 
in this temperature range are more strongly interacting than average, i.e. that these 
particles are part of agglomerates showing properties similar to those of the 17 vol% 
sample. 

4. Conclusion 

Extensive studies of the magnetic dynamics of a nanoparticle system containing nearly 
monodisperse ferromagnetic particles have been presented. We have shown that a 
strongly interacting particle system (the 17 vol% sample) displays critical dynamics 
reminiscent of that of a spin glass. Furthermore, the effect of an Arrhenius-Neel 
temperature dependence of r* has been explored and it has been found that dynamic 
scaling for this sample prevails although with slightly reduced values of the critical 
exponents zv and (5. 

For weakly interacting particle systems, the dipolar interaction will only slightly 
modify the relaxation compared to a non-interacting system, and the relaxation time 
can be obtained by introducing thermodynamical averages of the dipolar field in an 
expression for the single-particle relaxation time, as suggested in p3f . Such a theory 
will apply for the most dilute sample studied here (0.06 vol% sample) near the blocking 
temperature. 

For a wide range of particle concentrations (interaction strengths), neither a pure 
model for weakly interacting systems nor a model only assuming critical dynamics will 
correctly describe the magnetic relaxation of an interacting particle system. In the 
present work, this has been illustrated by the 5 vol% sample, in which correlations and 
collective behavior are of importance, as evidenced by non-equilibrium dynamics similar 
to that exhibited by spin glasses 0] and an apparent critical slowing down in a limited 
temperature range. Yet, the deviation from critical slowing down at lower temperatures 
and the failure to satisfy other signatures of critical dynamics show that the slowing 
down of the magnetic dynamics in this sample differs profoundly from that in spin 
glasses. The range of concentrations for which such a complex behaviour occurs becomes 
wider with increasing width of the particle size distribution, as can be evidenced by 
comparing to the system studied in JTTJ] . For the sample used in that study, the relative 
dipole interaction strength (= Mg(p/K, where is the volume fraction of particles) is 
comparable to that of the 17 vol% sample studied here, but that sample has a much 
wider energy barrier distribution and does not exhibit critical dynamics although it 
shows collective behaviour. 
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Figure 1. AC susceptibilities for the 0.06 vol% (□), 5 vol% (o), and 17 vol% (A) 
samples at frequencies of / =125 Hz (filled symbols) and / = 1000Hz (open symbols). 
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Figure 2. AC susceptibilities for the 5 vol% sample using the frequencies (left to right, 
open symbols) / =0.010, 0.031, 0.091, 0.31, 0.91, 3.1, 9.1, 31, 91, 310, 910, 3100, 9100 
Hz. The filled points are obtained from ZFC relaxation measurements and correspond 
to the frequencies / =0.019, 0.103, 0.979, 10.1 mHz. The dashed line indicates the 
equilibrium susceptibility. 
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Figure 3. AC susceptibilities for the 17 vol% sample using the frequencies (left to 
right) / =0.017, 0.051, 0.17, 0.51, 1.7, 5.1, 17, 51, 170 Hz. The dashed line indicates 
the equilibrium susceptibility. 




Figure 4. Relaxation rate for the 5 vol% sample obtained from ZFC relaxation 
measurements at the indicated temperatures after wait times of 300 s (open symbols) 
and 3000 s (filled symbols). 
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Figure 5. a) Relaxation time r c = ui^ 1 vs. reduced temperature for the 5 and 17 vol% 
samples. The lines are fits to equation (jij) as described in the text, b) Relaxation time 
vs temperature for the 0.06, 5 and 17 vol% samples. The freezing temperatures for the 
5 and 17 vol% samples were obtained from x' as described in the text. For the 0.06 vol% 
sample, the data correspond to the peak temperatures of \" \ an d the Arrhenius-Neel 
expression with parameter values as given in the text is shown as a dashed line. 
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Figure 6. Full dynamic scaling plot for the 17 vol% sample according to equation (^) 
using the Arrhenius-Neel expression for r* as described in the text. The asymptotic 
behaviours of G(x) are shown as straight lines. 
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Figure 7. (a) Plot of the AC data (open points) and ZFC relaxation data (filled 
points) according to equation ((?]) for the 5 vol% sample, (b) Plot of the AC data 
according to equation (0) for the 17 vol% sample. The dashed lines are obtained from 
G{x) with the expected asymptotic behaviours. Errorbars are only shown when larger 
than the symbol size. 




Figure 8. Relaxation rate vs. observation time at different temperatures. For the 
0.06 vol% sample (□) the relaxation rate is calculated from the AC susceptibility data. 
For the 5 vol% sample (o), the relaxation rate is calculated from AC susceptibility data 
(open symbols), ZFC data (filled symbols) and noise data (lines) and for the 17 vol% 
sample the relaxation rate is calculated from AC susceptibility data (A) and from G{x) 
with the expected asymptotic behaviours (dashed lines). 



